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Abstract 

i-^ ■ We investigate the effects of primordial non-Gaussianities in the primordial Uni- 

verse and the effects on the baryonic structure formation process. By relating 
the cosmic star formation rate in Gaussian and non-Gaussian scenarios to the de- 
tectability of high-redshift sources of reionization, we derive the expected Gamma- 
Ray Burst rate in the different models. We find that counts of high-redshift Gamma- 
Ray Bursts can be used as cosmological probes of non-Gaussianities and that they 
are suitable candidates to distinguish non-Gaussian effects at early epochs. 

1 Introduction 

Primordial density fluctuations grow during cosmological times and lead to the ob- 
served structures in the Universe [e.g. 1], via gas collapse and condensation in the 
evolving dark-matter haloes. The statistical distribution of the primordial fluctua- 
tions is still debated, though. In fact, despite the common assumption of Gaussian- 
ity, recent analyses show evidences for positively skewed distributions - see the 
summary Table 2 in [7]. 

High redshift represents an interesting regime to potentially investigate these non- 
Gaussian effects, since very early structures and primordial mini-haloes (hosting 
the first bursts of star formation, spreading the first heavy elements in the Universe, 
seeding the first black holes, being responsible for primordial cosmic re-ionization 
and hydrogen 21 -cm emission) are expected to be somehow affected by the under- 
lying matter distribution [8, 9, 16]. 

In this respect [8, 12], a key tool for studies of high-redshift environments might be 
represented by 7-ray bursts (GRBs), powerful explosions emitting 7 rays. These 
bursts have (comprehensive reviews are given in [17, 14]): isotropic equivalent 
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peak luminosities as high as ~ 10 erg s ; an isotropic angular distribution; and 
a bimodal duration distribution, made mostly of long GRBs (with a period longer 
than 2 seconds), and in minor part of short GRBs (detected at low redshift with 
a period shorter than 2 seconds). In the next, we will only consider long GRBs 
(LGRBs), because they are linked to the death of massive stars and hence are indi- 
cators of the local star formation episodes [e.g. 15, 2, 13]. 

Here, we will show that GRBs can be used as cosmological probe of the amount of 
non-Gaussianity present in the primordial density field, because they are sensitive 
to the underlying cosmological model through the cosmic star formation history - 
further datails on this subject can be found in [12]. 

We will reach this conclusion by performing a detailed analysis of the GRB rate 
in different Gaussian and non-Gaussian scenarios, based on N-body, hydrody- 
namic, chemistry simulations of early structure formation [7]. The star formation 
rate (SFR) and the adopted initial mass function (IMF) for the stellar populations 
tracked during the runs are used to determine the expected GRB formation rate den- 
sity in the various cases, and hence the integrated GRB rate (R), for both metal- 
poor population III (hereafter, popIII) regime and metal-enriched population II-I 
(hereafter, popII-I) regime. 



2 Model and results 

Given the cosmic SFR, p+, tracked by the different runs as a function of z in dif- 
ferent cosmological models [6, 7], the expected redshift distribution of observed 
GRBs can be computed once the GRB luminosity function (LF) and the GRB for- 
mation history have been specified [e.g. 19, 20, 21]. From these, we then compute 
the evolution of the GRB formation rate density, ti-grb? and the corresponding in- 
tegrated GRB rate, R. 

The (comoving) SFR densities are obtained by a set of numerical N-body, hydrody- 
namical, chemistry simulations with two different box sizes (0.5 and 100 Mpc//z, 
with h = 0.7) starting from initial conditions with a different level of primor- 
dial non-Gaussianity, parameterized in terms of the non-linear parameter /nl= 
0,10,50,100,1000. 

The observed peak photon flux, P, emitted by an isotropically radiating source at 
redshift z and corresponding luminosity distance cIl(z), as detected in the energy 
band E min < E < £ max , is 




(1) 
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where S(E) is the differential rest-frame photon luminosity of the source. To de- 
scribe the typical burst spectrum we adopt a "Band" function with low- and high- 
energy spectral indices equal to —1 and —2.25, respectively, and normalization: 

/•lOMeV 

L= E S(E) dE, (2) 

J IkeV 

with L isotropic-equivalent peak luminosity. 

Given a normalized GRB LF, ip(L), the observed number rate of bursts (in yr _1 ) 
at redshift z, with peak photon flux, P, between Pi and Pi is 



dt Jz dz' (l + z>) Jl Pi (z>) 



where the factor (1 + z') -1 accounts for cosmological time dilation, V(z') is the 
comoving volume at redshift z', and ncRB is the actual comoving GRB formation 
rate density. 

By assuming that GRBs are good tracers of star formation, we can write that the 
GRB formation rate density is proportional to the SFR density: 

«grb(z) = fcp*(z)> ( 4 ) 

where the normalization constant, k (whose dimensions are the inverse of a mass), 
incorporates further not-well-known effects, like GRB beaming, efficiencies, and 
black-hole production probability (depending on the adopted IMF). 
As for the LF, we will adopt a normalized GRB LF described by a single power-law 
with slope £ and decreasing exponentially below a characteristic luminosity, L+, 

Additionally, we take into account possible rdshift evolution of the GRB LF by 
setting = £*,o(l + Z Y '» where L* o is the characteristic luminosity at z = 0, 

and 5 is the evolution parameter. 

For simplicity, the normalization of ip(L) is included in k, and it is fixed when the 
GRB number rate in equation (3) is normalized to the rate observed at z = 0. 
From the previous relations we can finally compute the GRB rate (in units of 
yr _1 sr~ x ), R, as the derivative with respect to the solid angle, Q, of the GRB 
number rate, N: 
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Figure 1: Top panels: The expected popII-I and popIII GRB rates, R, in the 0.5 Mpc/h- 
side boxes (left), and in the 100 Mpc//i-side boxes (right), for models with different pri- 
mordial non-Gaussianity: /nl= (solid black line), 10 (dotted cyan line), 50 (short-dashed 
green line), 100 (long-dashed red line), and 1000 (dot-dashed blue line). Bottom panels: 
The corresponding relative contributions of the popIII GRB rates, Rm, to the total rate, R, 
for the different cosmological models. See [12]. 



We determine the parameters of the model by fitting the recent Swift redshift- 
complete sample by [21]. 

We apply our calculations to the different non-Gaussian cosmologies previously 
mentioned to derive count predictions as a function of z. We underline that the 
following results are valid for an ideal instrument, that is able to detect all the GRBs 
at any redshift. The main effects of primordial non-Gaussianities on GRBs are due 
to the differences in the redshift evolution of the SFRs in the various scenarios [see 
7] and, hence, in the consequently different GRB rates. 

In Fig. 1 (see also [12]), we plot the GRB rate, R, for the small 0.5 Mpc//j-size 
boxes (left panels) and the large 100 Mpc//i-size boxes (right panels). In the top 
panels, we show the redshift evolution for all the /nl values considered, while in 
the bottom panels we focus on the relative contribution of the popIII GRB rate 
(i?m) to the total rate, that is widely dominated by popII-I stellar generations, at 
redshift lower than z ~ 20. 

In both small and large volumes the effects due to the presence of primordial non- 
Gaussianities are visible at z > 10 — 15. The rates eventually converge at later 
times, when feedback processes start dominating the gas behaviour and the result- 
ing star formation. 
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These trends are valid for both popII-I and popIII regimes, even though the latter is 
usually negligible, predicting popIII GRB rates, Rm, that, following the behaviour 
of the popIII SFRs, drop by two orders of magnitude (bottom-left panel). Further 
discussions can be found in [12]. 

3 Conclusion 

Assuming that long 7-ray bursts are fair tracers of star formation [15, 2, 13], we 
propose to use them as probes of the underlying matter distribution at high red- 
shift, when the possible presence of non-Gaussianity would have visible effects on 
the baryon evolution (we refer the interested reader to [12] for further details and 
discussions). 

We find that already at z > 6 cosmologies with large /nl values have distinc- 
tive characteristics compared to those of the Gaussian case. At very early times 
(z ~ 15 — 20) the boost in the rate due to non-Gaussianities is ~ 2 — 3 orders of 
magnitudes for /nl=1000, and up to a factor of ~ 10 for /nl=100, while milder 
differences of a factor of ~ 2 are still visible for values of /nl~ 50. These effects 
are particularly evident on the total GRB rate, that is largely dominated by popII-I 
stars, while the result for the popIII GRB rate is noisier, mostly for /nl~ — 100, 
as a consequence of the corresponding, short-lived, popIII star forming regime 
[6, 10]. We highlight that differences in the GRB rate induced by non-Gaussianities 
are expected to be important mostly at very high redshift, while additional hydro- 
dynamic processes and feedback mechanisms would wash out any memory of pri- 
mordial non-Gaussianities at lower z. 

The existence of GRBs at such high redshift is not unlikely, as a small population 
of extremely dark GRBs, i.e. bursts for which the afterglow remains undetected in 
spite of early and deep Kband observations, has been recently identified [3]. While 
the nature of these GRBs is still matter of debate, it is possible that some of them 
are at z > 18. If confirmed, this could provide evidence for mildly positive non- 
Gaussian parameter. To draw more definitive conclusions and give more stringent 
constraints much larger high-z GRB complete samples, currently not available in 
the literature, are required. 

Acknowledgements 

This work has been performed under the HPC-EUROPA2 project (project num- 
ber: 228398) with the support of the European Commission - Capacities Area 
- Research Infrastructures. We also acknowledge the Italian Computing center 
(CINECA) for kind hospitality. 



5 



References 



. 9 
[20 



Barkana R., Loeb A., 2001, Phys. Rep., 349, 125 

Campisi M. A., Maio U., Salvaterra R., Ciardi B., 2011, MNRAS, 416, 2760 
D'Elia V., Stratta G., 2011, A&A, 532, A48 

Maio U., Dolag K., Meneghetti M., Moscardini L., Yoshida N., Baccigalupi 
C, Bartelmann M., Perrotta F., 2006, MNRAS, 373, 869 

Maio U., Dolag K., Ciardi B., Tomatore L., 2007, MNRAS, 379, 963 

Maio U., Ciardi B., Dolag K., Tornatore L., Khochfar S., 2010, MNRAS, 
407, 1003 

Maio U., Iannuzzi E, 2011, MNRAS, 415, 3021 

Maio U., 2011, Classical and Quantum Gravity, 28, 225015 

Maio U., Khochfar S., 2012, MNRAS, 421, 1113 

Maio U., Khochfar S., Johnson J. L., Ciardi B., 2011, MNRAS, 414, 1145 

Maio U., Koopmans L. V. E., Ciardi B., 2011, MNRAS, pp L197+ 

Maio, U., Salvaterra, R., Moscardini, L., & Ciardi, B. 2012, arXiv: 1209.2250 

Mannucci E, Salvaterra R., Campisi M. A., 2011, MNRAS, 414, 1263 

Meszaros P., 2006, Reports on Progress in Physics, 69, 2259 

Nuza S. E., Tissera P. B., Pellizza L. J., Lambas D. G., Scannapieco C, de 
Rossi M. E., 2007, MNRAS, 375, 665 

Petkova, M., & Maio, U. 2012, MNRAS, 422, 3067 

Piran T., 2004, Reviews of Modern Physics, 76, 1 143 

Salopek D. S., Bond J. R., 1990, Phys. Rev. D, 42, 3936 

Salvaterra R., Chincarini G., 2007, ApJ, 656, L49 

Salvaterra R., Delia Valle M., Campana S., Chincarini G, Covino S., 
D'Avanzo P., Fernandez-Soto A., Guidorzi C, et al 2009, Nature, 461, 1258 



[21] Salvaterra R., Campana S., Vergani S. D., Covino S., D'Avanzo P., Fugazza 
D., Ghirlanda G., Ghisellini G., Melandri A., Nava L., Sbarufatti B., Flores 
H., Piranomonte S., Tagliaferri G., 2012, ApJ, 749, 68 



6 



